Phase diagram and validity of one-parameter scaling near the two-dimensional 

metal-insulator transition 
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We explore the scaling description for a two-dimensional metal-insulator transition (MIT) of elec- 
trons in silicon. Near the MIT, Pt/p = (—l/p)d(lng)/dQnT) is universal (with p, a sample depen- 
dent exponent, determined separately; ^-conductance, T-temperature). We obtain the character- 
istic temperatures To and T\ demarking respectively the quantum critical region and the regime of 
validity of single parameter scaling in the metallic phase, and show that Ti vanishes as the transition 
is approached. For T < T\, the scaling of the data requires a second parameter. Moreover, all of 
the data can be described with two-parameter scaling at all densities - even far from the transition. 



PACS Nos. 71.30.+h, 73.40.Qv 

Many recent experiments have provided strong evi- 
dence for a two-dimensional (2D) disordered metallic 
state in Si metal-oxide-semiconductor field-effect tran- 
sistors (MOSFETs) and in other 2D systems 0-0, but 
the microscopic nature of the metallic phase still awaits 
a solid understanding. Since there have been predictions 
for 2D metallic behavior for quite some time j|] , it should 
not be too surprising to see a metal-insulator transition 
(MIT), but on the other hand it appears to violate rather 
general and widely accepted assumptions Hi. It is ap- 
parent from the experimental data PJ^JTp|-p^| that spin 
interactions are an important component of the metallic 
behavior and it is suspected that the dynam- 

ical Coulomb interactions are important j|,|l^|l4| since 
the non-interacting (screened) carriers are known to be 
insulating in 2D ||. 

The physics of the MIT can be reduced to certain scal- 
ing forms that describe the behavior of the conductivity 
<7 as a function of temperature T, sample length L, and 
carrier density n. A simple single-parameter scaling has 
proved to be a successful description of the critical be- 
havior of a in most 2D systems. It fails, however, in the 
metallic phase at the lowest temperatures, where u(T) 
becomes very weak and, in some systems, it even ap- 
pears to saturate as T — > 0. From the careful analysis 
of the experiments discussed below, we find that two- 
parameter scaling is more appropriate in that regime. 
In fact, we show that all of the data on both insulating 
and metallic sides of the transition can be described with 
two-parameter scaling at all densities - even far from the 
transition. Such scaling has been proposed theoretically 
near the MIT in the presence of dangerously irrelevant 
variables and it also provides an excellent descrip- 
tion of the MIT in a different type of the 2D metal [jL6|. 
This form of the two-parameter scaling function is consis- 
tent with the apparent success of the conventional single- 
parameter scaling in the limited range of (n,T) phase 
space. 

One signature of the 2D MIT is a reversal of the sign 



of da/dT. (a is written throughout dimcnsionlcssly - 
so all numerical values are multiplied implicitly by e 2 /h 
to recover SI units). Measurements of the conductance 
G = (w/L)a (w is device width) as a function of the 
carrier density n at several values of T reveal a crossing 
point at a critical density n c . Near enough to the MIT, 
a one-parameter scaling scheme holds, 



<r(n,T)~/(£) = /(<Wn> 



(1) 



where 5 n = (n — n c )/n c , <5* = S n /T p , and p — 1/zv is 
the scaling exponent describing the critical behavior of 
correlations scales of temporal and spatial fluctuations. 
f(8 n /TP) ~ exp^/TP) for T > T ~ \S n \ 1/p in the 
quantum critical region Jl3[ . Even for T < Tq, the data 
still can be scaled according to Eq. ([!]) to form a single 
curve. The single-parameter scaling fails at T = T\ < To, 
and we show that the scale T\ — Ti(S n ) vanishes as the 
transition is approached. 

We examine the phase diagram of a near the 2D MIT 
in two samples with peak mobility /Z4.2 — 0.5 m 2 /V sec 
at T = 4.2 K. To obtain T < 4.2 K, we used a He 3 - 
He 4 dilution refrigerator with heavily shielded wiring. 
Small AC signals were measured with lock- in amplifiers. 
The samples were generic two-probe Si MOSFETs with 
a 50 nm thick oxide layer; sample 17 was L x w = 5 x 
11.5 fim and sample 27 was 254 x 254 fim. Both devices 
have a density of oxide interface states N ox < 10 14 /m 2 , 
which is measured with standard techniques []l7f . 

Figure |l|(a) contains a(d n ,T) from a 2D inversion 
layer (sample 17) below 4 K, and the conventional one- 
parameter scaling [|l3) results appear in Fig. [|(b). The 
individual curves of a(d n ) at fixed T all cross at a crit- 
ical density n c as expected. As shown in Fig. |l|(b) the 
raw data for <r(<5*) collapse on to a single curve ["master 
curve" <7i(<$*) ~ f(Sn)} f° r a given 2D system if a certain 
exponent p is chosen in <5* = S n /T p . This collapse fits 
the suggestion of a quantum critical point at T = and 
5 n = 0. Critical conductivities in other systems fall in the 
range 0.3 < a c < 5, and so do the present results. Sample 
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FIG. 1. (a) ff(5 n , T) for 0.4 K < T < 4 K from sample 17 
with a schematic phase diagram inset. Arrows in the main fig- 
ure indicate the direction of flow for decreasing temperature. 
The phase diagram comprises two curves for To oc |(5, l | 1 / p 
(solid line) and Ti oc |c5 n | 1/ ' 9 (dotted line) that separate the 
S„ — T phase space into three regions. The quantum critical 
regime (QCR) is indicated schematically by cross-hatching. 
For T > Ti a single parameter can be used to scale the con- 
ductivity and for T < Ti scaling with one parameter fails, 
(b) Scaling of data in In a vs S n /T p for both samples (upper 
curve is sample 27 and lower curve is sample 17). The solid 
line emphasizes the linearity of In a in S n near the MIT. (c) 
In the regime of validity of one-parameter scaling (T > Ti), 
the universal behavior of zv/3t ~ vf3(g) is shown for six Si 
MOSFETs that span a wide range of parameters. 



27 exhibits the same features as sample 17 with slightly 
different parameters: a c of the two samples (0.32 for sam- 
ple 17 and 0.50 for sample 27) differ from each other by 
50% even though the values of n c (1.7 x 10 15 /m 2 and 
1.9 x 10 15 /m 2 , respectively) and the exponents p (1/2.0 
and 1/2.2, respectively) are within 10%. As predicted 
Q, the scaled data lna(S*) in Fig. 0(b) are approxi- 
mately linear near enough to the transition, for T > Tq 
(emphasized by the solid line). This is tantamount to 
the symmetry between a(5 n ) and l/a(—S n ), which has 
been predicted [jl3| and observed Q . Beyond the linear 
range (marked in the figure as To), the data continue to 
scale on both the insulating and metallic sides. At some 
point, however, the scaling fails on the metallic side and 
the individual traces of c(<5*) at fixed T diverge from 
<7i(<5*). These points in 5* can be read as positions in 
the phase diagram for the MIT. A phase diagram like 
that proposed for quantum phase transitions |]l9| is in- 
set in Fig. 0(a). The critical point is at T = and 
5 n = and crossover lines separate the quantum criti- 
cal region from the insulating and metallic regions. We 
have added a second (dotted) curve Tt(6 n ) demarking the 
region (T < Ti) where the conventional one-parameter 
scaling breaks down. T\ was obtained from individual 
traces of cr(<5*) by locating their departure from ai(£*) 
as shown in Fig. ||(a). It is obvious that the difference 
between one-parameter scaling and the experiment can 
reach as much as [cr(<5*) - ai(5*)]/a c = 35%. The T\ 
are plotted as a function of S n in Fig. ^(b), which shows 
clearly that the higher temperature data can be scaled 
over a wider range of 8 n . From Fig. ^(b), the character- 
istic temperature T± appears to have a power-law depen- 
dence on 8 n , as might be expected near a critical point. 
The power-law dependence of Ti obtains in both sam- 
ples and the exponents are very similar: 0.54 ± 0.01 for 
sample 17 and 0.49 ± 0.03 for sample 27. 

For T > Ti, one-parameter scaling holds. We show 
that, even though different parameters p, n c and a c are 
required to obtain the master curves, there is a univer- 
sal feature to the ultimate scaling function of the con- 
ductance as a function of length. In particular, we find 
that the function fir/p = (l/p)[— d(lng)/d(\nT)] is uni- 
versal. The results of this transformation are shown in 
Fig. ^(c) for the two samples studied here and compared 
with four other samples studied in the past [p|f20[| . Clearly 
the curve (3t/p is universal in the vicinity of the MIT. 
These experiments span two orders of magnitude in sam- 
ple dimension pQ] , an order of magnitude in critical ex- 
ponent p pof , about an order of magnitude in p.4.2, differ- 
ent surface electrostatics Q , and different manufacturing 
procedures. We expect the relation to hold in other 2D 
systems, and in fact a somewhat similar curve was ob- 
tained for higher peak-mobility Si MOSFETs With 
the assumption that the length over which quantum co- 
herence obtains is growing algebraically as T — > so 
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FIG. 2. (a) Method of observing crossover from 

one-parameter to two-parameter regimes by measuring de- 
viations from the o\ at each T. Dotted lines depict (arbitrar- 
ily chosen) 2% deviation, (b) Experimental phase diagram 
for both samples. Open and solid dots mark To, open and 
solid diamonds mark T\ . For each sample, the phase diagram 
is separated into the quantum critical regime (QCR) where 
lnu a <5*, Ti < T < To regime where the T-dependence of a 
is weaker but the one-parameter scaling still works, and the 
T < Ti regime where scaling can be achieved only with two 
parameters. 



/3(g) = d{hxg)/d{hxL) ~ dQjxg)/d(lnL (p ) as (3(g) ~ zp T . 
Fig. |](c) thus shows the universality of zv/3t ~ v (3(g). 
Recent numerical work has suggested Q that vfl(g) is 
a universal function close to the MIT. 

On the metallic side, for T < Ti, we find that each 
<r(<5*) trace deviates from the master curve in the same 
way, regardless of temperature, which suggests that it 
should be possible to scale all of the curves using two 
parameters. This, together with the predictions of two- 
parameter scaling for certain forms of MIT (e.g. with a 
strong triplet coupling, i.e. the spin-dependent part of 
the electron-electron interaction is large) JlfJ , encouraged 
us to examine a new scaling scheme according to 



(j(n,T) ~ T p f(S n /T p ). 



(2) 



This scaling form has been applied to a different form 
of the 2D metal with considerable success jl6|. As can 
be seen in Fig. ||, with this form we can scale all of 
the data for 0.4 < T < 1.5 K from the insulating side 
through to the metallic side using the same values of p 
that were used to accomplish single-parameter scaling at 
T > T\. For each sample, the data collapse onto a sin- 



gle (sample dependent) curve over a range of concentra- 
tions 6 n far larger than possible with the one-parameter 
scheme above. At T = 0.4 K, k B T w 0.02E F (E F - 
Fermi energy) at the MIT, i. e. lower than in many 
other MIT experiments on non-silicon systems. In addi- 
tion, Fig. U shows that the data obtained at even lower T 
(0.06 < T < 0.4 K) follow the same scaling law (Eq. |). 
For 5 n < 0.2, there is a maximum in o~(T) at T = T m and 
the scaling fails. This maximum has been attributed |l2] ] 
to Kondo coupling to disorder-induced local magnetic 
moments. The local moments did not affect the trans- 
port at T m < T, consistent with our finding (Fig. |^) 
that all of the data for T m < T follow the same two- 
parameter scaling form (Eq. ^|). Furthermore, this same 
two-parameter scaling holds even at T > 1.5 K if we al- 
low a small variation (< 5%) of n c with temperature (see 
Fig. H inset). This small change in n c might result from 
rescaled screening or other finite temperature effects . 

The power-law dependence of the prefactor for a on 
T is shown in the inset of Fig. [}| The new exponent p' 
appears with values 0.1 for sample 17 and 0.2 for sam- 
ple 27. Clearly, such a weak temperature dependence of 
the prefactor can become observable only when the scal- 
ing function f(S* ) itself becomes a very weak function of 
T, which happens at some T w Ti < To in the metallic 
regime. Indeed, Fig. |l|(c) (zv(3t) shows that, for aja c > 
4 in the metallic phase, /(<5*) becomes a power-law func- 
tion of T with an exponent comparable to the values of 
p' found from the two-parameter scaling. This explains 
naturally the apparent success of the single-parameter 
scaling in the description of the 2D MIT, and the exis- 
tence of the energy scale Ti(S n ) where it fails. The two 
samples, which differ in length by one and a half orders of 
magnitude, follow the scaling form (^|) albeit with some- 
what different exponents p' . It has been suggested |2(J 
that finite size effects may play a role in the failure of 
one-parameter scaling. On the other hand, a larger value 
of p' , i.e. a larger deviation from one-parameter scaling, 
is found in a larger sample, but it would be interesting 
to study if there is any systematic dependence of all the 
critical exponents on the sample length. 

The data (Fig. ||) also show that a c depends very 
weakly on temperature: a c — a(S n — 0, T) cx T p , re- 
sulting in <j c — at T = 0. It is important to note, how- 
ever, that the direct observation of a significant reduction 
of a c as T — > can not be achieved at experimentally 
accessible temperatures because of the very small expo- 
nent p' . For this reason, our result is not inconsistent 
with the constant values of o~ c observed in other experi- 
ments. Of course, we also can not distinguish this tem- 
perature dependence from a logarithm, which is another 
form proposed for corrections to scaling for certain quan- 
tum critical points J2l| . On the other hand, our result 
is in agreement with the recent work |l6| ] on the tran- 
sition between a different type of the 2D metal and an 
insulator, but the values of p' are considerably different, 
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FIG. 3. Two-parameter scaling of all of the data for both 
samples according to Eq. (2). The parameters from the new 
scaling law are inset. The critical concentration n c is com- 
pared to its low-temperature value to show the small change 
at higher temperature. For T < OAK only limited ranges of 
8„ can be scaled (see text). The arrows mark the end of the 
scaling curves if all data below 0.4 K are removed. 



presumably reflecting the different universality classes of 
the two situations. 

In the metallic phase, as T — > 0, a(T) becomes a very 
weak function after an initial rapid increase. As a result 
of this "saturation" ^] of o~(T), single-parameter scal- 
ing fails at the lowest temperatures (see also Ref. 28 ) 
and the scale T\ must exist. We have identified explicitly 
T\ = Ti(S n ) and shown that it vanishes as the transi- 
tion is approached. Furthermore, we have shown that all 
of the conductivity data in both metallic and insulating 
regimes can be scaled with two parameters according to 
Eq. (Q). That scaling form is consistent with the obser- 
vation that er(T) becomes very weak ("saturates") in the 
metallic phase as T —> 0. We reiterate that the scaling 
form (||) has been predicted for certain models . The 
predictions do not, however, decribe our data in detail 
(e.g. the sign or magnitude of p'). Scaling form (j|) has 
been also used successfully to describe a different type of 
a 2D MIT Jl|]. 
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